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THEORETICAL STUDY OF ORBITAL ELEMENT CONTROL WITH 


POTENTIAL APPLICATIONS TO MANNED SPACE MISSIONS 

By James C. Howard 

Ames Research Center 
Moffett Field, Calif. 

SUMMARY 


To facilitate the modification of established lunar or planetary orbits, a 
simple method has been developed for determining the impulsive velocity incre- 
ments required to make prescribed changes in the orbital elements. The direct 
determination of the velocity increments required to make prescribed changes in 
the elements of planetary orbits proceeds from a knowledge of the time rates of 
change of the orbital elements in the presence of perturbing forces. The sim- 
plification inherent in the present method involves the replacement of a set of 
differential equations relating the time rates of change of the orbital ele- 
ments to the components of the perturbing force by a set of ordinary equations 
relating the orbital element increments to the components of the impulse vector. 
These equations are then solved by conventional methods to obtain the impulse 
vector required to make prescribed changes in the orbital elements. This pro- 
cedure necessitates the omission of second and higher order quantities. How- 
ever, a 7090 digital computer program has been used to compute the errors 
involved, and the results obtained indicate that the errors are well within 
acceptable limits for the type of maneuvers contemplated. The equations have 
been mechanized on the 7 O 9 O digital computer and used to compute the impulsive 
velocity increments required to make independent changes in the orbital ele- 
ments of a variety of earth and lunar orbits. In addition, the velocity incre- 
ments required to produce orbital element changes of the type required to 
execute lunar reconnaissance maneuvers have been computed and the results pre- 
sented in graphical and tabular form. Specifically, the velocity increments 
required to modify a lunar orbit with a semimajor axis of 1,^00 statute miles 
and an eccentricity of 0.2 have been computed. Such an orbit would be of value 
if the object of a reconnaissance mission were to observe the lunar surface in 
the vicinity of perilune. Because the use of high eccentricity orbits can 
greatly reduce the velocity increments required to produce prescribed changes 
in orbit plane orientation, a lunar orbit with a semimajor axis of 2,200 stat- 
ute miles and an eccentricity of 0-5 has been used to further illustrate the 
method. It is assumed that such orbits would be established by successive 
applications of the method described in this report. 


INTRODUCTION 


Considerable research by Ames Research Center and other research groups 
has been directed towards reducing the complexity of guidance and control 



systems for manned space missions. References 1 and 2, for example, present 
results which describe some relatively simple navigation and guidance systems 
that have application to the midcourse return phase and the earth reentry phase 
of a manned lunar mission. As part of a continuing effort in this direction, 
an analytical study was undertaken to determine the possibility of devising 
simple methods of computing the impulsive-thrust-vector requirements for pro- 
ducing prescribed changes in the elements of an orbiting space vehicle. In 
order to determine the impulses required to change the size and shape of estab- 
lished orbits by the methods described in this report, it is necessary to know 
the semimajor axis, the eccentricity, and the true anomaly in the established 
orbit. In the case of earth orbits, the necessary orbital element data can be 
obtained from ground tracking stations. However, in the case of extraterres- 
trial missions, other back-up methods may be desirable. One method by which 
the semimajor axis and the eccentricity can be obtained is from a series of 
range determinations by the use of the parallax method. The determination of 
range by this method requires a knowledge of the radius of the target body and 
the angle subtended at the vehicle by this radius . Before true anomaly can be 
determined it is necessary to establish the position and direction of the line 
of apsides. However, discussion of the various methods of determining the 
orbit from an orbiting space vehicle is considered outside the scope of the 
present report, where it is assumed that the required orbital element data are 
available. This study was motivated, • in part, by the potential application of 
simplified computational procedures to space missions where fairly precise con- 
trol of the orbital elements is desired, for example, lunar or planetary recon- 
naissance missions, the shaping of parking orbits prior to injection into lunar 
or interplanetary trajectories, and the removal of the secular changes produced 
by the perturbations of the planet's gravitational field (ref. 3 ). The present 
report has two main objectives: 

(1) To describe a simple computational procedure for determining the 
impulsive velocity increments required to produce prescribed 
changes in the orbital elements of a space vehicle 

(2) To illustrate the potential application of this procedure to a 
lunar reconnaissance mission where prescribed orbital element 
changes may be made with relatively simple on-board computations. 

SYMBOLS 


semimajor axis, statute miles 

matrix element in the ith row and the jth column 
matrix element in the ith row and the jth column 
orbit eccentricity 




normal component of force 
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radial component of force 
transverse component of force 
universal constant of gravitation 
angular momentum per unit mass 

unit vector in the positive direction of the angular momentum vec- 
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orbit plane inclination 

a triad of inertially fixed, mutually orthogonal unit vectors 

(fig- 20) 

components of an impulse vector 

impulse required to make a prescribed change in the semimajor axis 


u ae 


impulse required to make prescribed changes in the semimajor axis 
and the eccentricity 
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impulse required to make a prescribed change in the eccentricity 

normal component of the impulse vector 

normal impulse required to make a prescribed change in the orbit 
plane inclination 

radial component of the impulse vector 

transverse component of the impulse vector 

radial impulse required to make prescribed changes in the argument 
of perifocus 

radial impulse required to null the changes induced in the argument 
of perifocus by I a 

radial impulse required to null the changes induced in the argument 
of perifocus by I ae 

radial impulse required to null the changes induced in the argument 
of perifocus by I e 

normal impulse required to make a prescribed change in the longitude 
of the ascending node 

semilatus rectum, a(l - e 2 ) 
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mass of space vehicle 

mass of the central body 

perturbing force vector 

position vector of vehicle in its orbit 

r 

unit vector in the positive direction of the position vector, — 
thrust vector 
true anomaly 

thrust direction measured from the local tangent to the flight path 


Pa 

direction 

of 

la 

P ae 

direction 

of 

^•ae 

Pe 

direction 

of 

Ie 

?a 

direction 

of 

■'-a 

7 ae 

direction 

of 

-*-ae 

7 e 

direction 

of 
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increment of the associated variable 

changes in the argument of perifocus associated with prescribed changes 
in the semimajor axis 

changes in the argument of perifocus associated with prescribed changes 
in the semimajor axis and the eccentricity 

changes in the argument of perifocus associated with prescribed changes 
in the orbit eccentricity 

2X2 determinant 

argument of latitude, to + v 

unit vector in the direction of increasing anomaly 
dynamical constant of gravitation, G(M + m) 
extraterrestrial gravitational potential function 
gravitational potential function for an oblate planet 
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argument of perifocus as measured from the ascending node line 
argument of perifocus, co + fi 
longitude of ascending node line 

unit vector in the direction of the ascending node line 
gradient of the gravitational potential function 
degrees of arc 

inverse of a matrix 
transpose of a matrix 


Subscripts 


components of subscripted impulse 

orbital element changes resulting in changes in the subscripted 
variable 


orbital element to be changed by the subscripted impulse 


represents the effects of extraterrestrial gravitation of subscripted 
potential function 

row and column, respectively, to which subscripted matrix element 
belongs 

represents the effects of oblateness of subscripted potential 
function 

radial, transverse, and normal components, respectively, of the 
subscripted variable 

orbital element to be changed by the impulse subscripted by a 

orbital element to be changed by the impulse subscripted by ae 

orbital element to be changed by the impulse subscripted by e 
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ANALYSIS 


The influence of thrust and gravitational forces on the elements that 
determine the size, shape, and orientation of the orbit of a space vehicle may- 
be determined by the use of Lagrange's planetary equations. The five elements 
that determine the size, shape, and orientation of a satellite orbit have to be 
controlled by a force vector having three components. It follows that the 
application of a force vector at an arbitrary point on an orbit will, in gen- 
eral, produce changes in all the orbital elements. For example, the applica- 
tion of a force vector which has a component normal to the orbit plane has the 
effect of producing a gyroscopic precession of the orbit plane; and changes in 
the orbit plane orientation give rise to changes in the argument of perifocus. 
Moreover, a force which lies wholly in the orbit plane has no influence on the 
orbit plane orientation, but does produce changes in the semimajor axis, the 
orbit eccentricity, and the argument of perifocus. If the differential equa- 
tions for the time rates of change of the orbital elements of a space vehicle 
in a non-Keplerian force field are considered, it is possible to make a direct 
determination of the impulses required to produce prescribed changes in the 
elements if second-order terns are neglected. If the second-order quantities 
are omitted the differential equations relating the time rates of change of the 
orbital elements to the components of the thrust vector can be replaced by a 
set of ordinary equations relating the orbital element increments to the com- 
ponents of the impulse vector. In addition to the influence of thrust which 
may be treated as a controlled perturbation, the equations for the time rates 
of change of the orbital elements can be used to study the effect of the earth's 
oblateness and extraterrestrial gravitational perturbations . If required, the 
terms representing the effect of the earth's oblateness can be used to deter- 
mine the changes produced in the elements of a parking orbit by the earth's 
gravitational perturbations. These changes, which depend on the length of time 
spent in orbit and on the orbit inclination, consist of periodic and secular 
variations. In considering the effect of the gravitational forces produced by 
an oblate planet having a potential function given by an infinite series of 
zonal harmonics, it is found that the greatest influence is exerted by the sec- 
ond harmonic of the planet's gravitational field. Although the influence of 
higher gravitational harmonics may be relatively smaller than the second, the 
cumulative effect of these harmonics is significant if the time spent in orbit 
is long enough to permit the secular changes to accumulate. 


Differential Equations 


A set of differential equations which may be used to study the influence 
of a perturbing force on the orbital elements of an orbiting space vehicle is 
derived in reference 4. See also references 5 and 6. For an alternative vec- 
torial derivation, see appendix B. These equations give the time rates of 
change of the orbital elements as functions of the radial, transverse, and the 
normal components of the perturbing force vector. Since the present study is 
concerned only with those elements that determine the size, shape, and orienta- 
tion of an orbit, the influence of perturbing forces on the time of perifocal 
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passage is omitted. The appropriate equations as derived in reference 4 are 
reproduced in appendix A. To conform to the most widely accepted space-age 
usage, changes of notation are introduced where necessary. The rates of change 
of the orbital elements are all expressed as functions of the true anomaly with 
the exception of the equation for the rate of change of the eccentricity. This 
equation is seen to depend on the true anomaly and the eccentric anomaly. For 
present applications, it is more convenient to have all the orbital element 
rates expressed as functions of the true anomaly only. For this reason, the 
eccentric anomaly is removed from the equation by expressing it as a function 
of the true anomaly and the eccentricity. The differential equations for the 
rates of change of the appropriate elements may be written in abbreviated form 
as follows : 


where 


Se ± 

e,- = 

i dF, J J 

ti 


i = 1, 2 ... 5 
j = 1 ; 2 , 3 


a " 1 ' a a ij( e i^ 




( 1 ) 


The repeated subscript on the right-hand side of equation (l) denotes a summa- 
tion. When equation (l) is written in full, the coefficients apj become the 
elements of a 5 X 3 matrix; eq are the five orbital elements that determine 
the size, shape, and orientation of the orbit and they are defined as follows: 
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The three mutually orthogonal components of the force vector are Fj. If sub- 
scripts r, Q, and h denote the radial, transverse, and normal components, 
respectively, the components Fi are defined by the following column vector: 
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( 3 ) 


^ 1 


1 F r 

f 2 

= 


Fa J 


l F h 


Equations For Impulsive Velocity Increments 


Assuming that second- order quantities may be neglected, the relationship 
between the orbital element increments and the components of impulse is given 

by 


Se i = dt 5t = a ij (F j 5t ) = a ±3 x 3 


(4) 


In terms of radial, transverse, and normal components, the impulse vector is 
defined as follows: 


11 ' 


1 Z r ' 

Is 
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I 3 j 


l ih 


In matrix notation, equation (4) assumes the following form: 


5a ^ 


j a ll a 12 0 ^ 

i 

5e 


a 21 a 22 0 


1 F r 

6to 

= 

a 31 a 32 a 33 


H 

CD 

Si 

i 


0 0 a 4 „ 

1 1 


l x h 

dQ j 


| 0 0 a53 j 



(5) 


( 6 ) 


It is clear from equation (6) that the problem of determining impulse consists 
in finding three unknown impulse components when there are five equations. 

Since it has five rows and three columns, the rank of the matrix which operates 
on the column vector of impulse components in equation (6) cannot exceed three. 
Hence, there cannot be more than three independent rows. It follows that any 
combination of these equations taken three at a time can be used to obtain an 
impulse vector having the magnitude and direction required to produce pre- 
scribed changes in the three corresponding elements. The influence of this 
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impulse on the remaining two orbital elements is obtained from the condition of 
compatibility or consistency. The arrangement of the zero elements in the 
matrix operator gives an indication of the way in which the orbital elements 
are influenced by the components of impulse. For example, the fact that 
a i 3 = a 23 = 0 denotes that changes in the semimajor axis and the eccentricity 
require in-plane forces. These two orbital elements are not influenced by 
forces normal to the orbit plane. Furthermore, the orbit plane inclination and 
the longitude of the ascending node line can only be changed by the application 
of forces normal to the orbit plane. This is evident from the fact that 


a 41 — a 42 - asi — as 2 — 0 

Since a 33 ^0, it follows that the argument of perifocus undergoes c han ges 
when either in-plane or out- of- plane forces are applied. If the argument of 
perifocus is defined as the sum of the nodal longitude and the argument of the 
latitude of perifocus, then changes in orbit plane orientation give rise to 
changes in the argument of perifocus- The matrix element a 33 is a measure of 
the change in the argument of perifocus associated with changes in the orbit 
plane orientation. The matrix elements aqj are functions of the true anomaly 
and the argument of the latitude; hence, some of these elements vanish at cer- 
tain orbital locations . The orbital locations where a given matrix element 
vanishes are points of ineffectiveness for the corresponding impulse component. 
However, such locations can be used to advantage if it is required to manipu- 
late one orbital element without changing the others. To show the dependence 
of the matrix elements on the true anomaly and the argument of the latitude, 
equation (6) is written in full. 



( 7 ) 
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The matrix element a 43 which controls the orbit plane inclination is seen to 
vanish when the argument of the latitude is 90° or 270° , whereas a 53 , which 
controls the longitude of the ascending node line,, vanishes at the ascending 
and descending nodes. Hence , a normal impulse applied at a nodal position will 
alter the orbit plane inclination without changing the longitude of the ascend- 
ing node line; whereas, a normal impulse applied when the argument of the lati- 
tude is 90° or 270° will alter the longitude of the ascending node line without 
changing the orbit plane inclination. The components of the impulse vector 
required to make prescribed changes in the semimajor axis, the eccentricity, and 
the argument of perifocus can be determined from the following matrix equation: 


Therefore 


a n a i 2 0 ' 


1 I r | 


1 &a 

8.21 ^22 0 



= 

Be 

a 31 a 32 a 33 / 


1 Ih ) 


1 


M 


1 a n 

a 12 

0 ' 

-1 

1 ^ 

I© 

= 

a 21 

a 22 

0 


&e 

x h J 


^ a 31 

a 3 2 

a 33 / 


1 


The impulse vector obtained from equation (9) would produce the prescribed 
changes in the orbital elements; however, such an impulse would produce 
unwanted changes in the orbit plane orientation. Subsequent corrections to the 
orbit plane orientation would then re-introduce errors in the argument of peri- 
focus. In view of this it is considered more expedient to effect the pre- 
scribed changes by the use of an impulse vector which lies wholly in the orbit 
plane. The influence of such an impulse on the semimajor axis, the eccentric- 
ity, and the argument of perifocus can be determined from the following matrix 
equation: 


a ll a 12 1 

1 1 

1 



Hr ! 



a 21 a 22 



= 

6e 



HeJ 


1 

, a 3l a 32 / 


\ 


( 10 ) 


The rank of the matrix which operates on the column vector of impulse compo- 
nents cannot exceed two, and therefore the number of linearly independent rows 
cannot exceed two. Hence, any combination of these equations taken two at a 
time can be used to obtain an impulse having the magnitude and direction 
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required to produce prescribed changes in the corresponding elements. The 
impulse vector required to produce prescribed changes in the semimajor axis and 
the eccentricity is obtained from equation (10) in the following form: 


' x r \ 

( a n 

a 12 ) 

- 1 / 5a ' 

1*9 1 1 

( a 21 

a 22 j 

! \ be J 


( 11 ) 


Compatibility or consistency of equation (lO) requires that 


Therefore , 




( a 31 a 32 ) 


l 


f a n 

a 12 

\ a 21 

a 22 


-i 


' Sa 
( be 


a ll 

a 12 

a 21 

a 22 

a 31 

3-32 



/ 


5a 

be 



a 32 ) 


^ a ll 

a 12 \ 1 

f 8a ] 

a 21 

a 22 | 

Se / 


( 12 ) 


(13) 


Physically, this means that in producing the prescribed changes 6a and be by 

the application of the impulse vector (l r Ig) , the change given by equation 
(12) has been induced in the argument of perifocus. From equation (7) it is 
seen that these induced errors can be removed without changing the remaining 
elements if a radial impulse is applied at perifocus or apofocus. In the event 
that the impulse vector should have an out-of -plane component Ift, equation (8) 
shows that the total change induced in the argument of perifocus is given by 


5w = 


a 31 a 32 


a ll 

a 12 \ “I 

! 5a | 

a 21 

a 22 / 

\ be j 


+ a 33 Ih 


( 1*0 


The solution of equation (ll) 


is given by 
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1 _ 1 

/ a 22 

“ a 12 

5a 


\ - a 2i 

a ll j ^ 

^ be 

A = (an 

a 22 - 

a 12 a 21 ) 



( 15 ) 


(16) 


On substitution from equation (7) , it is found that the impulse vector required 
to produce prescribed changes 5a and 5e in the semimajor axis and the eccen- 
tricity, respectively, is given by 



(17) 


Radial impulse requ ired to change the arg ument of perifocus . - If the radial 
impulse required to produce an incremental change Sto in the argument of peri- 
focus be denoted by I w , then from equation (7) 



If this impulse were applied at perifocus or apofocus, it would have no effect 
on the remaining orbital elements. Hence, the radial impulse required to 
change the argument of perifocus by 6w without changing the remaining orbital 
elements is given by 

^ S co (l8a) 

Sco 

where 



(l8b) 
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Normal Impulses required to change the orientation of the orbital plane. - 
If the impulse required to produce an incremental change Si in the orbit 
plane inclination be denoted by Ii, then from equation ( 7 ) 


li = 


M- (l + e cos v) 

i(l - e 2 ) cos(v + u) 


Si 


(19) 


Likewise, if the impulse required to produce an incremental change S ft in the 
longitude of the ascending node be denoted by Iq, then from equation ( 7 ) 


= 


P (l + e cos v) 
a(l - e 2 ) sin(v + w) 


sin 1 




( 20 ) 


From equations ( 19 ) and (20), it is seen that if a normal impulse were applied 
when the argument of the latitude (v + co) was 0° or l 80 °, the orbit plane 
inclination would be changed without altering the longitude of the ascending 
node. Alternatively, if it were required to alter the longitude of the ascend- 
ing node without changing the orbit plane inclination, the normal impulse would 
be applied when the argument of the latitude was 90 ° or 270 °. 


Impulsive Velocity Increments Required to Produce Prescribed 
Changes in the Orbital Elements of Noncircular Orbits 


Simultaneous cha nges in the semimajor axis and the eccentricity.- To 
facilitate the computation of the velocity increments required to produce pre- 
scribed changes in the semimajor axis and the eccentricity of the orbit of a 
space vehicle, equation ( 15 ) may be rearranged as follows: 


/M 


1 1 

j aa 22 

-eai2 \ 


1 — 

' a 

CD 

H 


A 

1 ~ 9.9-21 

eau J 


^ Be 


( 21 ) 


The elements of this matrix may be obtained from equation ( 17 ). By letting 
5a/a = rn/lOO and Se/e = n/100, we can use equation (21) to obtain the impulse 
vector required to produce an m percentage change in the semimajor axis, and 
an n percentage change in the orbit eccentricity. The required impulse vec- 
tor has the following components: 


1 Ir 


1 Cn + C 12 \ 

CD 

H 


C 2 1 + c 22 j 


( 22 ) 
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where the matrix elements 


Cij have the following values. 


• 


" 

m 

l ^ 

cos v(2 + e cos v) + e | 

200 

N 

! a(l - e 2 ) 

sin v 1 


( 23 ) 


'12 - 


ne 

/ 11 

fd + e cos v) 2 1 

100 

N 

1 a(l - e 2 ) 3 

sin v _J 


(24) 


C21 


m / l-i ~ 
200 y a (l - e2) 


(l + e cos v) 


(25) 


O22 


ne 2 / I- 1 

1°° V a(l - e 2 ) 3 


(l + e 


cos v) 


(26) 


The magnitude of the impulse vector is given by 

lae = J (O11 + C 12 ) 2 + (C 21 + C 22 ) 2 


To achieve the desired result, the impulse vector must be inclined at an angle 

7 to the local horizontal, where 
/ ae ’ 


7 


ae 


tan 


On + 012^ 
C 2 i + C 22J 


(27) 


The angle between the impulse vector and the local tangent to the flight path 
at the point of application of the impulse is given by P ae , where 


3 ae 


= tan 1 


/Cn + C 12 \ 
\C 2 i + C 22 y 


tan 


e sin v 
1 + e cos v 


(28) 
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Changes in the argument of perifocus resulting from changes in the semi- 
major axis a nd the orbit eccentricity .- If the change in the argument of peri- 
focus induced by the application of the impulse Iae be denoted by 6co ae , then 


& w ae - ( a 3i a 32) 



(29) 


where 


a 31 


1 / a(l - e a T 

e 7 P 


(30) 


and 


a 32 


1 / a(l - e a ) sin v( 2 •+ e cos v) 

e J ^ (l + e cos v) 


(31) 


A radial impulse is required to null S w without changing the remaining 

orbital elements, and it must be applied at perifocus or apofocus. If this 

impulse be denoted by I w , then from equation (l8a) 

8.0 


"^ae 


aico 

<3w 


S u ae 


(32) 


Prescribed cha nges in the semimajor axis . - To produce an m percentage 
change in the semimajor axis of a noncircular orbit, without changing the orbit 
eccentricity, requires the application of an impulse I a , where 


la = jCn 2 + C 21 2 


(33) 


For the desired effect, the impulse vector must be applied at an angle 7 a to 
the local horizontal where 


7 


a 



(34) 
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The angle between the impulse vector and the local tangent to the flight path 
at the point of application of the impulse is given by 3 , where 

cl 


Pa = 


tan 




sin v 
e cos v 


(35) 


Changes in the argu ment o f p erifocus resulting from changes in the semi- 
major axis . - If the change in the argument of perifocus induced by the impulse 
I be denoted by 5 to a , then 


S u a = ( a 3i 8-32) ( J (36) 

w 


where a 31 and a 32 have the values given in equations (3O) and (31). As indi- 
cated previously, a radial impulse is required to null Sw a without changing 
the remaining orbital elements and it must be applied at perifocus or apofocus. 
If the required impulse be denoted by 1 ^ , then from equation (l8a), 

Si 



dim 

doo 


Sw a 


(37) 


Prescribed changes in t he orbit eccentricity. - The impulse required to 
produce an n percentage change in the eccentricity, without changing the 
semimajor axis, is given by 


(38) 

The impulse vector must be applied at an angle / e to the local horizontal 
where 


I* = 


-I 


'12 


+ c 


22 


7 e 



(39) 
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The corresponding angle between the impulse vector and the local tangent to the 
flight path is given by 


Pe 



- tan 


-l 


e sin v 
1 + e cos v 


(40) 


From equations (24) , ( 26 ), and (40) it is seen that 


C lg f e sin v A 
C 22 \1 + e cos v/ 


(41) 


hence , the angle (3 is always ±90°. It follows that to change the orbit 
eccentricity without changing the semimajor axis, the impulse must be directed 
along the normal to the flight path. 


Changes in the argument of perifocus resulting f rom changes in the orbit 
eccentricity . - The change in the argument of perifocus produced by the impulse 
I e is given by 


f Gi A 

6 We = ( a 3i a 32) ( ] (42) 

\C 22 J 


The radial impulse required to null Swe is given by 



5lw 

da) 


6w e 


(43) 


Application to Lunar Reconnaissance Maneuvers 


Since only a limited area of the lunar surface can be observed from one 
orbit, the orbit plane orientation must be changed to permit observation of 
other areas. Orbital plane rotation about any position vector can be accom- 
plished by the application of an impulse normal to the orbital plane at any 
desired orbital location, as indicated in previous sections. 

The orientation of an orbital plane is by convention defined in terms of 
the orbit plane inclination and the longitude of the ascending node line. This 
way of resolving the orbit plane orientation into components is convenient for 
establishing the orientation of the plane in space. However, for lunar or 
planetary reconnaissance work, there is no need to determine the orbital plane 
orientation in terms of the angular displacement of the node line from some 
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fixed line, and the angular displacement of the orbital plane from some fixed 
plane, measured as a rotation about the node line. For the applications con- 
templated here, it is more convenient to consider the orbital plane rotation 
vector rather than its components. It is assumed that from some initial orbit 
which is unsatisfactory from the point of view of observation and photographic 
survey, an orbit with a satisfactory size and shape has been established by the 
methods previously indicated. When the surface of the target body has been 
observed and photographed from this orbit, the orbital plane can then be 
rotated about the position vector from the center of the target body to the 
vehicle by the application of an impulse normal to the orbital plane. This is 
evident from the following considerations: the preces_sional torque required to 

change the direction of the angular momentum vector h without changing its 
magnitude is given by 


M = a) X h 


(WO 


where M is the precessional torque vector and w is the precessional rate 
vector. This symbol is not to be confused with the symbol for the argument of 
perifocus. Solving this equation for to gives 


GJ 


h x M 
h 2 


however, 


M = r X T 


(45) 

(46) 


where r is the position vector and T is the thrust vector. On substituting 
from equation (46) in equation (45), the angular velocity cj is obtained in 
the following form: 


- h X (r X T) F h r 

O) = i - 


h c 


h 


A 

r 


(47) 


where F^ is the component of thrust normal to the orbit plane. Equation (47) 
shows that the orbital plane rotates about the position vector in the presence 
of a force which has a component in the direction of the normal to the orbital 
plane. It is seen that for a given force F^, the angular velocity u> is pro- 
portional to the radial distance r; and, hence, the impulse required to pro- 
duce a prescribed change in orientation decreases with radial distance. 

The normal impulse 1^ required to produce a prescribed change So in 
the orbit plane orientation can be obtained from equation (47) as follows: 
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CO 


Fhr 

h 


A 

r 


do A 
dt r 


therefore 


50 = i st = I ( F * 5t ) - e x h 


therefore 


Ih = 7 50 


I 


M- 


l( 1 - e 2 ) 


(l + e cos v) 


So 


(48) 


RESULTS AND DISCUSSION 


The results presented here have been divided into two parts: one part 

deals with the velocity increments required to make prescribed changes in the 
elements of earth orbits, and the other deals with the application of the 
method to lunar reconnaissance maneuvers. 


Earth Orbital Element Changes 


The results plotted in figures 1 and 2 show the velocity increments 
required to produce a range of prescribed changes in the semimajor axis of an 
earth orbit, subject to the constraint that the eccentricity retains its ini- 
tial value. The results plotted in figures 1 and 2 are for an orbit with a 
semimajor axis of 8,120 statute miles. However, the curves in figure 1 refer 
to an orbit with an eccentricity of 0.2; whereas those of figure 2 refer to an 
orbit with an eccentricity of 0-5 • Each figure shows the velocity increments 
required to produce five percentage increases in the semimajor axis. The per- 
centage increases selected for computation were: 0.1, 0.3, 0.5, 0 . 7 , and 0.9* 

For convenience, however, the corresponding increase in miles is noted on each 
curve. The curves superimposed on these plots show the directions which the 
velocity increments must have in order to produce the desired change. It is 
interesting to note that there are two points on each orbit where a tangential 
impulse can be used to produce a prescribed change in the semimajor axis with- 
out changing the eccentricity of the orbit. It is easy to show that this con- 
dition arises when the length of the position vector to the vehicle equals the 
length of the semimajor axis. The condition r = a is satisfied when the true 
anomaly satisfies the following equation: 


v 


cos -1 ( -e) 
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Changes in the a rgument of perigee resulting from prescribed changes in 
the semimajor axis.- The results plotted in figures 3 and 4 show the changes in 
the argument of perigee resulting from the range of percentage increases in the 
semimajor axis of the earth orbit discussed in the preceding section. Each 
curve of figure 3 shows the changes in the argument of perigee resulting from 
the corresponding curve of velocity increments plotted in figure 1. Likewise, 
each curve of figure 4 shows the changes in the argument of perigee resulting 
from the corresponding curve of velocity increments plotted in figure 2 . 

Eccentricity changes . - The impulsive velocity increments required to pro- 
duce prescribed changes in the orbit eccentricity subject to the constraint 
that the semimajor axis retains its initial value are plotted in figures 5 and 
6 . The results shown are for an earth orbit with a semimajor axis of 8,120 
miles. However, the curves shown in figure 5 refer to an orbit with an eccen- 
tricity of 0 . 2 ; whereas those of figure 6 refer to an orbit with an eccentric- 
ity of 0.5- Each figure shows the velocity increments required to produce five 
percentage increases in the orbit eccentricity. The percentage increases 
selected for computation were: 0.1, 0 . 3 , 0 . 5 , 0 . 7 , and 0.9* The curves super- 

imposed on these plots show the directions the velocity increments must have in 
order to produce the required change in the eccentricity. To change the eccen- 
tricity without changing the semimajor axis, it has been shown in equation (4l) 
that it is necessary to apply the impulse vector in a direction perpendicular 
to the local flight path, and in the plane of the orbit. See the curve of |3 
in figures 5 and 6 . 

Changes in the argume nt of perige e resulting from ^res_cribed_changes_ in 
the orbit eccentricity .- The curves plotted in figures 7 and 8" show the incre- 
mental changes induced in the argument of perigee by the velocity increments 
used to produce the percentage changes in eccentricity discussed in the pre- 
vious section. Each curve of figure 7 shows the changes in the argument of 
perigee induced by the corresponding curve of velocity increments plotted in 
figure 5 . Likewise, each curve of figure 8 shows the changes in the argument 
of perigee induced by the velocity increments plotted in figure 6. It is of 
interest to note that there are two points on an orbit where a normal impulse 
can be applied to produce a prescribed change in the eccentricity without induc- 
ing a change in the argument of perigee. It is easy to show that these two 
points correspond to values of the true anomaly given by the following equation: 


v = l 80 ° ± tan 



Orbital plane inclina t ion .- To produce a prescribed change in the orbital 
plane inclination, a velocity increment must be applied in the direction of the 
normal to the orbit plane. Equation (19) has been used to compute the magni- 
tude of the impulsive velocity increments required to produce a 1 ° change in 
the orbit plane inclination of an earth orbit. The results are plotted in fig- 
ure 9 . The initial values of the orbital elements are as indicated on the 
figure. 


20 



Nodal longitude . - As in the case of orbit plane inclination, changes in 
the nodal longitude require that the corrective impulses be applied in the 
direction of the normal to the orbital plane. Equation (20) has been used to 
compute the impulsive velocity increments required to produce a 1° change in 
the longitude of the ascending node line of an earth orbit with an inclination 
of 45° • The results are plotted in figure 10. The initial values of the 
orbital elements are as indicated on the figure. It is to be noted that 
changes in the orbit plane inclination are coupled to changes in the longitude 
of the ascending node line, and vice versa, unless the impulses are applied 
when the argument of the latitude (v + oj) is 0°, 90 °, l80°, or 270 °. 


Application to Lunar Reconnaissance 


Advantage of using high eccentricity orbits for certain types of recon - 
naissance maneuvers .- It is evident from equation (48) that fuel requirements 
could be reduced for certain types of reconnaissance missions if the lunar sur- 
face area of interest were in the vicinity of perilune and subsequently the 
orbit plane were rotated about the line of apsides. The impulse required to 
produce a l u change in orbit plane orientation is seen to be a minimum if the 
orbit plane is made to rotate about the line of apsides by the application of a 
normal impulse at apolune . The impulse required is a maximum if the plane 
change maneuver is executed at perilune. The ratio of the impulse required at 
apolune to that required at perilune is given by (l - e)/(l + e). Equation 
(48) has been used to compute the velocity increments required to rotate vari- 
ous orbit planes through a 1° angle about the line of apsides. The results in 
table I give a good indication of the advantage of executing plane change 
maneuvers when the vehicle is at apolune of high eccentricity orbits. 

Prescribed changes in the semimajor axis of a lunar orbit . - Equations (33) 
through (35) have been used to compute the impulse vector required to make pre- 
scribed changes in the semimajor axis of a lunar orbit, subject to the con- 
straint that Ae = 0. The results for a lunar orbit with a semimajor axis of 
1,600 miles are plotted in figures 11 and 12. The orbit eccentricities have 
the values indicated on the figures. The incremental changes in the argument 
of perilune resulting from changes in the semimajor axis are plotted in figures 
13 and 14. 

Changes in perilune height .- A lunar orbit with a semimajor axis of 1,400 
miles and an eccentricity of 0.2 has a perilune height of 40 miles. The veloc- 
ity increments required to reduce the perilune height by 5.6 miles were com- 
puted and plotted as a function of the true anomaly in figure 15* Curve 1 
shows the velocity increments required to produce the required change in height 
by reducing the semimajor axis without changing the orbit eccentricity. Curve 
1 of figure l6 shows the changes in the argument of perilune induced by changes 
in the semimajor axis. When the perilune height is changed without altering 
the orbital period (i.e., without changing the semimajor axis), curve 2 of fig- 
ure 15 results. The corresponding plot of induced changes in the argument of 
perilune is given by curve 2 of figure l6. Because the use of high eccentric- 
ity orbits can greatly reduce the velocity increments required to produce 
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prescribed changes in orbit plane orientation, the velocity increments required 
to modify a lunar orbit with a semimajor axis of 2,200 miles and an eccentric- 
ity of 0.5 have been computed. It is assumed that this orbit which has a peri- 
lune height of 20 miles would be established by successive applications of the 
method described in this report. The results plotted in figure 17 show the 
velocity increments required to reduce the height at perilune by 5>5 miles. 
Curve 1 of figure 17 shows the velocity increments required to reduce perilune 
height by reducing the semimajor axis without changing the orbit eccentricity. 
Curve 2 of figure 17 shows the velocity increments required to modify the 
height at perilune by the prescribed amount, by changing the orbit eccentricity 
without altering the semimajor axis. When equal percentage changes are made in 
the semimajor axis and the orbit eccentricity, curve 3 of figure 17 results. 

The corresponding curves of figure 18 show the changes induced in the argument 
of perilune by the velocity increments plotted in figure 17 - The curves of 
figure 19 give the directions which the velocity increments plotted in figure 
17 must have in order to produce the prescribed change in perilune height. 


CONCLUDING REMARKS 


Replacing a set of differential equations relating the time rates of 
change of the orbital elements to the components of a perturbing force with a 
set of ordinary equations relating the orbital element increments to the com- 
ponents of an impulse vector makes it possible to compute the impulsive veloc- 
ity increments required for prescribed changes in the elements of lunar or 
planetary orbits. The necessity of omitting second and higher order quantities, 
which the method entails, gives rise to errors in the computed results. How- 
ever, numerical analysis indicates that these errors are well within acceptable 
limits for the types of maneuvers contemplated. As an example, a lunar orbit 
with a semimajor axis of 1,600 miles and an eccentricity of 0.2 may be consid- 
ered. In making a 5 - mil e change in the semimajor axis, it was found that the 
impulse computed on the basis of the simplifying assumptions made gave rise to 
a change in the semimajor axis which was approximately 100 feet in error. As a 
second example, an earth orbit with a semimajor axis of 8,000 miles and an 
eccentricity of 0.1 may be cited. In making an 8-mile change in the semimajor 
axis, it was found that the computed impulse gave rise of an incremental change 
in the semimajor axis which was approximately 50 feet in error. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 3 , 196^ 
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APPENDIX A 


DIFFERENTIAL EQUATIONS FOR THE TIME RATES OF CHANGE OF THE 
ORBITAL ELEMENTS IN THE PRESENCE OF PERTURBING FORCES 


It is shown in reference ij- that the variation of the orbital elements can 
be expressed in terms of the components of the disturbing force. The formula- 
tion is such that the rates of change of the orbital elements are expressed as 
functions of the radial, transverse, and normal components of the perturbing 
force vector. In terms of the notation of reference k, the equations for the 
rates of change of the orbital elements assume the following form: 


da 

dt 



[S tan cp sin w + T sec cp(l + e cos w) ] 


(Al) 


de 

dt 



cos cp [S sin w + T(cos w + cos E) ] 


(A 2 ) 


^ q | -aS cos 2 9 cos w +rT sin w (2 + e cos w) +rW sincp tan[ (l/ 2 )i]sin u} 
dt - sin 9 cos 9 


(A3) 


di _ rW cos u 
dt Jpa cos 9 


(A 4 ) 


dG 

dt 


rW sin u 

Jp a cos 9 sin i 


(A5) 


where 

w true anomaly 

sin 9 eccentricity 
E eccentric anomaly 

u arg ume nt of the latitude (i.e., the sum of the argument of perifocus, 

measured from the ascending node line and the true anomaly) 

S,T,W radial, transverse, and normal components, respectively, of the per- 
turbing force vector 
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In the notation of the present report (see fig. 20) 

S = r • (Vcp ob + Vcp et + T) (A6) 

T = 0 • (Vcp ob + Vcp et + T) (A7) 

W = h • (Vcp Qb + Vcp et + T) (A 8 ) 

u = CO + v (A9) 

w = v (A10) 

sin 9 = e (All) 


The rates of change of the orbital elements are all expressed as functions of 
the true anomaly, with the exception of the equation for the rate of change of 
eccentricity. This equation is seen to depend on the true anomaly and the 
eccentric anomaly. For present applications it is more convenient to have all 
the orbital element rates expressed as functions of the true anomaly only. For 
this reason the eccentric anomaly is removed from the equation by expressing it 
as a function of the true anomaly and the eccentricity as follows: 


a(l - e 2 ) 

r = — 

1 + e cos v 



e cos E 


therefore 


cos E 


cos v + e 
1 + e cos v 


therefore 


cos v + cos E 


cos v(2 + e cos v) + e 
1 + e cos v 


(A12) 


2k 



Substituting from equations (A 6) through (A12) in equations (Al) through (A5) 
yields the following form of the equations for the rates of change of the 
orbital elements: 


da _ 2a 2 

^ «7|j.a(l - e 2 ) 


[(e 


sin v)r + (l + 


e 


cos v)$] * (v9 ob + V9 e t + T ) 


(A13) 


de /a(l - e 2 ) f / . \^ , 

at = 7 ^ — t (s “ T)r 


cos v(2 + e cos v) + e 
1 + e cos v 


e j- * (vq>ob + vq> et + T ) 

(All)-) 


dco _ /a(l - e 2 ) 

dt J q 


( 


1 9* , 

— cos v j r + 
e 


1 sin v(2 + e cos v) 
e 1 + e cos v 


A 

e 


“n ( V t f:! ( i /g)l! } «) • <Wob + Wet ♦ T) 


(A15) 


di 

dt 


/a(l - e 2 ) cos(v + to) A N 

J s L 1 + e W V h • Wob + Wet + T) 


(Al6) 


da 

dt 


/a(l - e 2 ) 

sin i(l + e cos v) 


sin(v + w) A 

1 l h 


(vcpob + V9et + T ) 


(A17) 


where 


„ Stp A ,1 9® ^ , 1 

Vcp = ^-*-r + — t- 1 -0 + 


Sep A 


dr r 9u r sin u di 


h 


(A18) 
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APPENDIX B 


VECTOR DERIVATION OF THE EQUATIONS FOR 
ORBITAL ELEMENT PERTURBATIONS 


The equation of motion of a particle of unit mass moving in an inverse - 
square -law central force field is 


d 2 r _ jj_ A 
dt 2 r 2 P 


(Bl) 


Vector multiplication of each side of equation (Bl) by the angular momentum 
vector h gives 



X h = - — ■ (r X h) 
r 


(B2) 


h = r x v = (r 2 0)h 


(B3) 


On substitution from equation (B3) in equation (B2) it is seen that 


d 2 r 

dt 2 


x h = (p0) 


A 

9 


therefore 


d v X h 
dt q 


• A d A 

99 = — r 
dt 


(B4) 


The integral of equation (B4) is given by 


v X h A _ 

■ = r + e 

M- 


(B5) 


where e is a constant vector of integration. The vector e may be expressed 
in terms of its scalar magnitude and a vector of unit length as follows: 


e = e'a (B 6) 

where a is a unit vector (fig. 20). Substituting equation (b 6) in equation 
(B5) gives 
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(B7) 


v x h 

P 


A , A 
r + ea 


Equation (B7) may be solved to obtain the position vector r. S£alar multi- 
plication of each side by "r gives the following equation for r: 


r • (v X h) 
P 


r + 



therefore 


(r X v) • h 


= r(l + e cos v) 


therefore 


r = 


h 2 /p 


1 + e cos v 


(B8) 


Vector multiplication of each side of equation (B7) by h gives 


h X (v X h) = ph X (ea + r) 


therefore 


v = ^ [0 + e(h X a) ] 


(B9) 


With the notation of figure 20, the unit vector a may be expressed in the 
following form: 


A , vA 

a = (cos v)r 


(sin v)0 


(BIO) 


If the assumption of an inverse-square-law central force field is not satis- 
fied, the equation of motion must be modified accordingly. In the presence of 
a perturbing force P, the equation of motion becomes 


(Bll) 


Furthermore, in the presence of the perturbing force vector P, the assumption 
of constancy no longer applies to the vector e. Hence, 


dv _ p P£ 

dt r 3 
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(B12) 




de 

dt 



x h + v x (r x P) 



[ (r X v) X r] 


therefore 


V- 


de 

dt 




X h + v X (r X P) 


therefore 


11 ~ = P X h + v X (rxP) 
dt 


(B13) 


The first term on the right side of equation (BI 3 ) may be written in the fol- 
lowing alternative form: 


P X h = h(re - 0r) • P 


(B14) 


likewise, 


— ,AA AA 

X P = r(h6 - 0h) 


P 


(B15) 


When substitutions are made from equations (B9) and (B15), the second term on 
the right side of equation (BI 3 ) is given by 


v X (r X P) 


_ qr 
h 


AA 

r 0 + e 


aA 

a 9 + 


,A 

(0 


A,AA 

a)hh 


• P 


(Bl6) 


Substitution for a from, equation (BIO), equation (Bl6) becomes 

v X (r X P) = ^~ [(l + e cos v)r0 - e sin v(e^ + hh) ] • P (B17) 
From equations (Bl4) and (B17) it follows that 


Pxh+vX (r X P) = 


(l + e cos v) + h 


aA 

re 


ure aa , Aa\ 

Tj— sin v 00 + her ) 


qre 

h 


aa9 

sin v hhr 


P 


(B18) 
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The components of the term on the left side of equation (B13) may be obtained 
as follows : 


M- 


de 

dt 



co X 



(B19) 


where co is the angular velocity of the vector e- A As the vector e" moves in 
the orbit plane , it rotates about the unit vector h which is normal to the 
plane. In addition, the orbit plane rotates about the instantaneous position 
vector. Hence, the vector co is given by 


therefore 


co 


coh + 


h x (r x P) 


h c 


co 


A 

coh + 


r 

h 


,aA 

(rh) 


P 


(B20) 


By substitution from equation (B20) in equation (B19) it is found that 


de 

dt 


, . A f . . \ A 

p [ e cos v + coe sm v J r + p I coe cos v - e sin v J 0 


pre 

h 


AA. 


sin v(hh) • P 


(B21) 


Equating coefficients of like vectors in equations (Bl8) and (B21) yields the 
following results: 


and 


p (^e cos v + coe sin v^ = 2h(§ • P) 


pre . a , a \ — 

p [ coe cos v - e sin v ) = - ( — sxn v0 + hr ) • P 


On solving for e and (!) it is found that 




cos v(2 + e cos v) + 


A ( h . \ Al 

9 + I — sm v r t • P 


(B22) 


co = 


r / 

t sin v(2 + e cos v)0 
eh ' 


h \ A 

— cos v r 
pe 


. P 


(B2b 
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It is seen that for circular orbits there is a singularity in the equation for 
(I) because of the presence of e in the denominator. The physical signifi- 
cance of this is that for low eccentricity orbits the argument of perifocus is 
not well defined and does not exist for circular orbits. 


Orbital Plane Orientation 


Since the orientation of a plane in space is uniquely determined by the 
normal to its surface, the orientation of an orbit plane is determined by the 
angular momentum vector. This fact may be used to advantage in finding the 
ti me rates of change of the elements defining the orientation of an orbit plane 
in space. The longitude of the ascending node line A and the inclination of the 
orbit plane are related to the unit vectors k and h as follows: 


A A , .A 

k X h = (sin i)ft 


therefore 


d^ 


di \ £ 


\ 


dt 


. dft 


kX If = (cos i £)B + (sin i (kX ft) 


\ 


dt 


A 


(B24) 


and 


A dh r /A A, /A — 

k X ^ | (k X 0) (h • P) 


(B25) 


A 


The unit vector 9 has the following components along and perpendicular to the 
ascending node line 

$ = -(sin 0$ + cos 9 (6 X ft) (B2 6) 


On substitution from 
obtained: 

A 

k x 

From equations (B24) 


equation (B26) in equation (B25), the following result is 

(h • P) (B27) 


A 

dh _ r 
dt h 


,A A, A 

sin 0(k X ft) + cos 6 cos ifl 


and ( B27 ) it follows that 


di r cos 9 /A — 

at " ~ — (h ‘ p) 


(B28) 


dft r sin 6 /A — . 

= Z — • (h • P) 

dt h sm i 


(B29) 
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Orbital Plane Orientation and the Argument of Perifocus 


Since the orbital plane rotates about the instantaneous position vector r, 
the argument of perifocus as measured from the line of nodes will vary in the 
manner indicated in sketch (a). 


Zk 



Sketch (a) 


When the sine rule is applied to the spherical triangle ABC, it is seen that 


sin 9 

sin(i + Bi) 


sin(e + 58) 
sin i 


therefore 


B0 = -(tan 0 cot i)5i 


Since the position vector is assumed constant during the orbital plane rota- 
tion, it follows that 6 9 = Boo. Hence, 


dco 

dt 


-(tan 9 cot i) 


di 

dt 


(B30) 
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Substituting for di/dt from equation (B 28 ), the rate of change of the argu- 
ment of perifocus produced by the rotation of the orbit plane is given by 


dco r cot i sin 9 A — . 

at h (h ' p) 


(B31) 


If the argument of perifocus is measured from the inertially fixed x axis, 
then the rate of change associated with orientation changes is given by the 
following sum: 


dfi dw 
dt dt 


(B 32 ) 


Substituting for dfi/dt from equation (B 29 ) in equation (B32), the required 
rate is obtained in the following form: 


dfi dm r sin Q tan[(l/2)i] A p-. 

dt dt h 1 ’ ' 


(B 32 a) 


If equations (B32a) and (B23) are combined, the total rate of change of the 
argument of perifocus is given by 


dm 

dt 


h 

pe 


COS V 


A 

r 


+ sin v(2 + 
eh 


e cos 


A 

v)e 


+ 


r sin (v + to) tan [(l/2) i] A 


h 


}• 


(B33) 


The semimajor axis is defined as follows: 

Let 

(h • h)/p 

d — ' 

1 - e • e 


(B34) 


It is known from orbit theory that for closed orbits, the quantity 2 a is the 
distance from perifocus to apofocus . The quantity a is the semimajor axis. 
In an ideal inverse -square -law central force field it remains constant. How- 
ever, the force field assumed in the present study is noncentral so that the 
quantity a must be treated as a variable. From equation (B 9 ), the square of 
velocity is given by 



therefore 

1 = 2 _ vf 

a r p 


(B 35 ) 


(B36) 
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Differentiating equation (B36) with respect to time gives the following rate 
of change of a: 


1 da 2 dr _ 2v dv 

a 2 dt r 2 dt p dt 


therefore 


_1_ da = _2_ f , a _ dr 4 2 (- _ dv 

a 2 dt r 2 \ dty p \ dt 


therefore 


a.' 


2 dt r‘ 


^a = A f* . dr ) + 2 - . fp _ Br 


dty p 


therefore 


da _ 2a /— — 

dt " 1 


(v • P) 


(B37) 


Substituting for v from equation (B 9 ) yields 


da _ 2a g 
dt " h 


[(e sin v)r + (l + e cos v)^] • P 


(B38) 


Perturbing Forces 


Perturbing forces are assumed to be generated by the earth's gravity field 
and extraterrestrial gravity fields. Both of these forces are assumed to be 
derivable from known gravitational potential functions. It is further assumed 
that thrust forces are available to introduce controlled perturbations. In 
terms of_the gravity gradients and the thrust forces, the perturbing force 
vector P may be expressed as follows: 


P = Vcp ob + W e t + T 


(B39) 


For convenience of reference the results are summarized below: 


da 2a c 


dt h 


[(e sin v)r + (l + e cos v)e] • (Vcp ob + Vcp e t + T) 


(B 40 ) 


de _ J / h . \ a i r 

dt |\P Sin V r h 


cos v(2 + e cos v) + 


9 f ’ (V 9 ob + Wet + T ) 


(B 41 ) 
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da3 

dt 


(- 


h | A r 

— cos v r + -T- 
He J eh 


sin v(2 + e cos v) 


A 

e 


r sin(v + co)tan[(l/2)i]l M , — , 

1 £ U J |h j • (vcp ob + V9 e t + T ) 


(B42) 


di r cos(v + u) A 
dt h n 


(vcpob + W e t + T ) 


(B43) 


d2 _ r sin(v + 10 ) A 
dt h sin i 


(Wob + W e t + T ) 


(B44) 
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TABLE I.- VELOCITY INCREMENTS REQUIRED TO ROTATE THE LUNAR-ORBIT PLANE THROUGH 


1° ABOUT THE LINE OF APSIDES 




Impulse, 

a, 


ft/sec 

miles 






Perifocus 

Apofocus 



Velocity increment, ft/sec 
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Figure 1.- Velocity increments required to make prescribed changes in the semimajor axis of an earth 

satellite subject to the constraint that Ae = 0. 
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Figure 2.- Velocity increments required to make prescribed changes in the semimajor axis of an earth 

satellite subject to the constraint that Ae = 0. 
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Figure 4.- Changes in the argument of perigee resulting from prescribed changes in the semimajor axis. 
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Figure 5 *- Velocity increments required to make prescribed changes in the eccentricity of an earth 

satellite subject to the constraint that Aa = 0. 
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Figure 6. - Velocity increments required to make prescribed changes in the eccentricity of an earth 

satellite subject to the constraint that Aa = 0. 
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Figure 7 • - Changes in the argument of perigee resulting from prescribed changes in the orbit 

eccentricity. 
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Figure 8 
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- Changes in the argument of perigee resulting from prescribed changed, in the orbit 

eccentricity. 
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Figure 9*- Velocity increments required to produce a 1° change in the orbit plane inclination of an 

earth orbit. 
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Figure 11.- Velocity increments required to make prescribed changes in the semimajor axis of a lunar 

orbit subject to the constraint that Ae = 0. 
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Figure 12.- Velocity increments required to make prescribed changes in the semimajor axis of a lunar 

orbit subject to the constraint that Ae = 0. 
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Figure l4.- Changes in the argument of perilune resulting from prescribed changes in the semimajor axis. 
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Figure 15- - Impulse vector required to reduce perilune height by 5.6 miles. 
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Figure 19- - Thrust direction for velocity increme 
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